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Intergalactic Magnetic Field and Arrival Direction of 
Ultra-High-Energy Protons 

Dongsu Ryu 1 ' 4 , Santabrata Das 2 , and Hyesung Kang 3 
ABSTRACT 

We studied how the intergalactic magnetic field (IGMF) affects the propaga- 
tion of super-GZK protons that originate from extragalactic sources within the 
local GZK sphere. Toward this end, we set up hypothetical sources of ultra-high- 
energy cosmic-rays (UHECRs), virtual observers, and the magnetized cosmic web 
in a model universe constructed from cosmological structure formation simula- 
tions. We then arranged a set of reference objects mimicking active galactic nuclei 
(AGNs) in the local universe, with which correlations of simulated UHECR events 
are analyzed. With our model IGMF, the deflection angle between the arrival 
direction of super-GZK protons and the sky position of their actual sources is 
quite large with the mean value of (9) ~ 15° and the median value of 9 ~ 7 — 10°. 
On the other hand, the separation angle between the arrival direction and the sky 
position of nearest reference objects is substantially smaller with (5*) ~ 3.5 — 4°, 
which is similar to the mean angular distance in the sky to nearest neighbors 
among the reference objects. This is a direct consequence of our model that the 
sources, observers, reference objects, and the IGMF all trace the matter distribu- 
tion of the universe. The result implies that extragalactic objects lying closest to 
the arrival direction of UHECRs are not necessary their actual sources. With our 
model for the distribution of reference objects, the fraction of super-GZK pro- 
ton events, whose closest AGNs are true sources, is less than 1/3. We discussed 
implications of our findings for correlation studies of real UHECR events. 
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Introduction 



The nature and origin of ultra-high-energy cosmic rays (UHECRs), especially above the 
so-called Greisen-Zatsepin-Kuz'min (GZK) energy of E GZ k ~ 50 EeV (1 EeV = 10 18 eV), 
has been one of mo st perplexing puzzles in a strophysics over five decades and still remains to 
be understood (see iNagano fc Watsonll2000l . for a review). The highes t energy CR dete cted 
so far is the Fly's Eye event with an estimated energy of ~ 300 EeV (jf3ird et al. 1994 ). At 
these high energies protons and nuclei cannot be confined and accelerated effectively within 
our Galaxy, so the sources of UHECRs are likely to be extragalactic. 

At energies higher than E GZ k, it is expected that protons lose energy and nuclei are 
photo-disintegrated via the interactions with the cosmi c microwave background radiation 
(CMB) along their traje ctories in the intergalactic space (jGreisenll 19661 ; IZatsepin &: Kuz'min 
19661 ; iPuget et al.lll976l ). The former is known as the GZK effect. So a significant suppres- 
sion in the energy spectrum above -Egzk could be regarded as an obse rvational evidence 
for the extragalactic origin of UHECRs (see, e.g., iBerezinsky et al.ll2006l ). However, the ac- 
curate measurement of the UHECR spectrum is very difficult, partly because of extremely 
low flux of UHECRs. But a more serious hurdle is the uncertainties in the energy calibra- 
tion i nherent in det ecting and modeling extensive airshower events (e.g.. INagano fc Watson 
2000l ; IWatsonll2006l ). Nevertheless, both the Yakutsk Extensive Air Shower Array (Yakutsk) 
and t he High Resolution Fly's Eyes (HiRes) reported observations of the GZK suppres- 
sion (jEgorova et al.l 120041 ; lAbbasi et al.l l2008al ) , while t he Akeno Giant Air Shower Array 
(AGASA) claimed a conflicting finding of no suppression (jShinozaki fc Teshimall2004l ). More 
recent data from the Pierre Auger Observatory (Auger) support the existence of the GZK 
suppression ([Abraham et al.ll2008bl : ISchussler et al.ll2009l ). Below E GZ k, however, the four 
experiments reported the fluxes that are different from each other by up to a facto r of several, 
imply ing the possible existence of systematic errors in their energy calibrations (IBerezinsky 
2009h . 



The overall sky distribution of the arriv al directions of UHECRs below Eqzk seems 
to support the isotropy hypothesis (see, e.g., INagano fc Watsonl bood ). This is consistent 
with the expectation of uniform distribution of extragalactic sources; the interaction length 
(i.e. horizon distance) of protons below Eqzk is a few Gpc and the universe can be consid- 
ered homogeneous and isotropic on such a large scale. The horizon distance for super-GZK 
events, however, decreases sharply with energy and -Rqzk ~ 100 Mpc for E = lOOEeV 
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( jBerezinsky fc Grigor'eval Il988l ). The matter distribution inside the local GZK horizon 
(Rqzk) is inhomogeneous. Since powerful astronomical objects are likely to form at deep 
gravitational potential wells, we expect the distribution of the UHECR sources would be 
inhomogeneous as well. Hence, if super-GZK proton events point their sources, their arrival 
directions should be anisotropic. 

The anisotropy of super-GZK events, hence, has been regarded to provide an important 
clue that unveils the sources of UHECRs. So far, however, the claims derived from analy- 
ses of different experiments are often tantalizing and sometimes conflicting. For instance, 
with an excessive number of pairs and one triplet in the arrival directi on of CRs above 



40 EeV, the AGASA dat 


1996 




Takeda et al. 


1999) 



Takeda et al.lll999l ). On t he other hand, the HiRes stereo data are consistent with 



the hypothesis of null clustering (jAbbasi et al.ll2004l . |2009| ). The auto-correlation analys is of 
the Auger data reported a weak excess of pairs for E > 57 EeV ([Abraham et al.ll2008al ). In 
addition, the Auger Collaboration found a correlation between highest energy events and the 
large scale structure (LSS) o f the uni verse using nearby active galactic nuclei (AGNs) in the 



2008a 



Veron-Cettv Veronl (20061 ) catalog (jThe Pierre Auger Collaborationll2007l : Abraham et al. 



_ Hague et al.l 120091 ) as well as using nearby objects in different catalogs (jAublin et al. 
20091 ) . A correlation between highest AGASA events with nearby galaxies from SDSS was 
reported (ITakami et al.l 120091 ). The HiRe s data, however, do not show such correlation of 
highest energy events with nearby AGNs (jAbbasi et al.ll2008bl ). but instead show a correla- 
tion with distant BL Lac objects (jAbbasi et al.l 120061 ). 



The interpretation of anisotropy and correlation analyses is, however, complicated owing 
to the intervening galactic magnetic field (GMF) and intergalactic magnetic field (IGMF); 
the trajectories of UHECRs are deflected by the magnetic fields as they propagate through 
the space between sources and us, and hence, their arrival directions are altered. Even with 
considerable observational and theoretical efforts, however, the nature of the GMF and the 
IGMF is still poorly constrained. Yet, models for the GMF generally assume a strength of 
~ a few At G and a coherence length of ~ 1 kpc for the field in the Galactic halo (see, e.g., 
Stanevlll997l ). an d predict the deflectio n of UHE protons due to the GMF to be 9 ~ a few 
degrees (see, e.g., ITakami k, Satol 120081 ) . The situation for the IGMF in the LSS has been 
confusing. Adopting a mo del for the IGMF with the average strength of (B) ~ 100 nG in 
filaments, ISigl et al.l (120031 ) showed that the deflection of UHECRs due to the IGMF could be 
very large, e.g., 9 > 20° for protons above 100 EeV. On the other hand, iDolag et al.l (120051 ) 
adopted a model with (B) ~ 0.1 nG in filaments and showed that the deflection should be 
negligible, e.g., 9 <C 1° for protons with 100 EeV. 



Recently, iRyu et al.l (j2008l ) proposed a physically motivated model for the IGMF, in 
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which a part of the gravitational energy released during structure formation is transferred to 
the magnetic field energy as a result of the turbulent dynamo amplification of weak seed fields 
in the LSS of the universe. In the model, the IGMF follows largely the matter distribution 
in the cosmic web, and the strength and coherence length are predicted to be (B) ~ 10 nG 
and ~ 1 Mpc for the field in filam ents. Such field in filaments is expect ed to induce th e 
Faraday rotation (ICho fc Rvul 120091 ) , whi ch is consistent with observation (IXu et al.l 120061 ). 



With this model IGMF, iDas et al.l (j2008l ) (Paper I hereafter) calculated the trajectories of 
UHE protons (E > 10 EeV) that were injected at extragalactic sources associated with 
the LSS in a simulated model universe. We then estimated that only ~ 35 % of UHE 
protons above 60 EeV would arrive at us with 8 < 5° and the average value of deflection 
angle would be (8) ~ 15°. Note that the deflection angle of (8) ~ 15° is much larger than 
the angular window of 3.1° used by the Auger collaboratio n in the study of the correlation 
betw e en highest energy UHECR events and ne arby AGNs (IThe Pierre Auger Collaboration 



20071 : Abraham et al.ll2008al : Hague et al.l 120091 ). 



In this contribution, as a follow-up work of Paper I, we investigate the effects of the 
IGMF on the arrival direction of super-GZK protons above 60 EeV coming from sources 
within 75 Mpc. The limiting parameters for energy and source distance are chosen to match 
the recent analysis of the Auger collaboration. Without knowing the true sources of UHE- 
CRs, the statistics that can be obtained with observational data from experiments are lim- 
ited; some statistics that are essential to reveal the nature of sources are difficult or even 
impossible to be constructed. On the other hand, with data from simulations, any statis- 
tics can be explore d. In that sense, simulations complement experiments. Here, with the 
IGMF suggested by lRyu et al.l (120081 ). we argue that the large deflection angle of super-GZK 
protons due to the IGMF is not inconsistent with the anisotropy and correlation recently 
reported by the Auger collaboration. However, the large deflection angle implies that the 
nearest object to a UHECR event in the sky is not necessarily its actual source. In Section 
2, we describe our models for the LSS of the universe, IGMF, observers, and sources of 
UHECRs, reference objects for correlation study, and simulations. In Section 3, we present 
the results, followed by a summary and discussion in Section 4. 



2. Models and Simulations 

In our study, the following elements are necessary: 1) a model for the IGMF on the LSS, 
2) a set of virtual observers that represent "us", an observer at the Earth, in a statistical 
way, 3) a set of hypothetical sources of UHE protons with a specified injection spectrum, 
and 4) a set of reference objects with which we performed a correlation study of simulated 
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events. In Paper I, we described in detail how we set up 1), 2), and 3) by using data from 
cosmological structure formation simulations. Below, we briefly summarize models for 1), 
2), and 3) and explain in details the reason to introduce "reference objects" in this study. 



2.1. Large Scale Structure of the Universe 

We assumed a concordance ACDM model with the following parameters: Obm = 0.043, 
Q BM = 0.227, and tt A = 0.73, h = F /(100 km/s/Mpc) = 0.7, and a 8 = 0.8. The model 
universe for the LSS was generated through simulations in a cubic region of comoving size 
100/i _1 (= 143) Mpc with 512 3 grid zones for gas and gravity and 256 3 particles f o r dar k 



matter, using a PM/Eulerian hydrodynamic cosmology code described in lRyu et al.l (119931 ). 
The simulations have a uniform spatial resolution of 195.3/i -1 kpc. The standard set of 
gasdynamic variables, the gas density, p g , temperature, T, and the flow velocity, v, were 
used to calculate the quantities required in our model such as the X-ray emission weighted 
temperature T x , the vorticity, u, and the turbulent energy density, e tU rb, at each grid. 



2.2. Intergalactic Magnetic Field 



We adopted the IGMF from the model by iRyu et al.l (120081 ): the model proposes that 
turbulent-flow motions are induced via the cascade of the vorticity generated at cosmological 
shocks during the formation of the LSS of the universe, and the IGMF is produced as a 
consequence of the amplification of weak seed fields of any origin by the turbulence. Then, 
the energy density (or the strength) of the IGMF can be estimated with the eddy turnover 
number and the turbulent energy density as follow: 



£ B = <P | ~ | ^turb- 

v teddy , 

Here, the eddy turnover time is defined as the reciprocal of the vorticity at driving scales, 
t eddy = I/deriving (oj = V x v) , and is the conversion factor from turbulent to magnetic 
energy that depends on the eddy turnover number t/t eddy . The eddy turnover number was 
estimated as the age of universe times the magnitude of the local vorticity, that is, t age u>. The 
local vorticity and turbulent energy density were calculated from cosmological simulations for 
structure formation described above. A functional form for the conversion factor was derived 
from a separate, incompressible, magnetohydrodynamic (MHD) simulation of turbulence 
dynamo. 



For the direction of the IGMF, we used that of the passive fields from cosmological 
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simul ations, in which magnetic fields were generated through the Biermann battery mecha- 
nism (lBiermann||l950) at cosmologica l shocks and evolved passively along with flow motions 
( IKulsrud et al.lll997l ; iRyu et al.lll998l ). In principle, if we had performed full MHD simula- 
tions, we could have followed the amplification of the IGMF through turbulence dynamo. In 
practice, however, the currently available computational resources do not allow a numerical 
resolution high enough to reproduce the full development of MHD turbulence. Since the 
numerical resistivity is larger than the physical resistivity by many orders of magnitude, 
the growt h of magnetic fields is saturated before dynamo action becomes full y operative 



(see, e .g., IKulsrud et al.l 1 199 71 ). This is the reason why we adopted the model of IRyu et al. 



( 120081 ) to estimate the strength of the IGMF, but we still used the the passive fields from 
cosmological simulations to model the field direction. 

Figure 1 shows the distribution of magnetic field strength in a slice of (143 Mpc) 2 in 
our model universe. It shows that the IGMF is structured like the matter in the cosmic 
web. As a matter of fact, the distribution of the IGMF is very well correlated with that 
of matter. The strongest magnetic field of B > 0.1/xG is found in and around clusters, 
while the field is weaker in filaments, sheets, and voids. In filaments which are mostly 
composed of the warm- hot intergalactic medium (WHIM) with T = 10 5 - 10 7 K, the IGMF 
has (B) ~ 10 nG and (B 2 } 1/2 ~ a few x 10 nG JRvu et al.l I2~0~0~8m . Note that the deflection 
of UHECRs arises mostly due to the field in filaments (see Paper I). The energy density of 
the IGMF in filaments is Eb ~ 10~ 16 ergs cm" 3 , which is a few times smaller than the gas 
thermal energy density and an order of magnitude smaller than the gas kinetic energy density 
thereQ The IGMF in filaments induces the Faraday rotation; the root-mean-squar e (rms) 



value of rotation measure (RM) is predicted to be a few rad m 1 (ICho fc Ryull2009l ). That 
is consisten t with th e value s of RM toward the Hercules and Perseus- Pisces superclusters 
reported in ku et all (bood lfl 



2.3. Observer Locations 

In the study of the arrival direction of UHECRs, the IGMF around us, that is, in the 
Local Group, is important too. It would have been ideal to place "the observer" where the 



1 We note that our mo del does not include a possible contribution to the IGMF from galactic black holes, 
AGN feedback (see, e.g.. iKronberg et al.ll2001l ); so our model may be regarded to provide a baseline for the 
IGMF. With the contribution, the real IGMF might be even stronger, resulting in even larger deflection (see 
Section 3.1). 



2 The values of |RM| i n Xu et al.l (|2006l ) is an order of magnitude larger than the value above, a few rad 



m . However. IXu et al.l ( 20061 ) quoted the path length, which is about two orders of magnitude larger. 
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IGMF is similar to that in the Local group. Unfortunately, however, little is known about the 
IGMF in the Local Group. Hence, instead, we placed "virtual" observers based on the X-ray 
emission weighted temperature Tx- The groups of galaxies t hat have the halo temperatur e 
similar to that of the Local Group, 0.05 keV < kT x < 0.5 keV (IRasmussen fc Pedersenll200ll ). 
were identified. About 1400 observer locations were chosen by the temperature criterion. In 
reality, there should be only one observer on the Earth. But in our modeling we could choose 
a number of observer locations to represent statistically "us" without loss of generality, since 
the simulated universe is only one statistical representation of the real universe. Then, we 
modeled observers as a sphere of radius 0.5/i" 1 Mpc located at the center of host groups, 
in order to reduce the computing time to a practical level. The distribution of handful 
observers are shown schematically in Figure 1. One can see that the observers (groups) are 
not distributed uniformly, but instead they are located mostly along filaments. 



2.4. AGNs as Reference Objects 



As noted in Introduction, the Auger Collaboration recently reported a correlation be- 
tween the direction of t heir highest energy events an d the sky position of nearby AGNs from 
the 12th edition of the IVeron-Cetty &: Veronl (120061 ) (VCV) catalog; the correlation has the 
maximum significance for UHECRs wi t h E > 60 EeV and AGNs with distance D < 75 Mpc 
(IThe Pierre Auger Collaboration!l2007l ; Abraham et al.ll2008al ; lHague et al.ll2009l ). There are 
about 450 AGNs with D < 75 Mpc (more precisely, 442 AGNs with redshift z < 0.018, for 
which the maximum significance of the correlation was found) in the VCV catalog. In that 
study it is not known which subclass of those AGNs or what fraction of them are really 
true sources of UHECRS. Here, we regard those AGNs as "reference objects", with which 
correlation studies are performed. 

In order to compare our correlation study with that of the Auger collaboration, we 
specified the following condition to determine "model" reference objects in the simulated 
universe: 1) the number of the objects within 75 Mpc from each observer should be on 
average ~ 500, and 2) their spatial distribution should trace the LSS in a way similar to the 
AGN distribution in the real universe. To set up the location of such reference objects, we 
identified "clusters" of galaxies with kTx ^0.1 keV in the simulated universe. Of course 
some of these clusters with kT x < a few keV should be classified as groups of galaxies. But 
for simplicity we call all of them as clusters. The reason behind this selection condition is 
that the gas temperature is directly related with the depth of gravitational potential well; 
the hottest gas resides in the densest, most nonlinear regions of the LSS where the most 
luminous and energetic objects (e.g. AGNs) form through frequent mergers of galaxies. We 
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then assumed that each cluster hosts one reference object at its center. For each observer, 
we generated a list of reference objects inside a sphere of radius 75 Mpc, whose number 
is on average ~ 500; the exact number of reference objects varies somewhat for different 
observers. Then, each observer has its own sky distribution of reference objects, with which 
we studied the correlation of simulated events. Although our reference objects could be any 
astronomical objects that trace the LSS, hereafter we refer them as "model AGNs", because 
the selection criteria were chosen to match the number of AGNs with that from the VCV 
catalog. Figure 1 shows the schematic distribution of handful model AGNs at the center 
of host clusters. Obviously the host clusters (and the AGNs) are not uniformly distributed 
either. 

In our set-up, the distance to the model AGNs, D, can be arbitrarily small. In reality, 
however, the closest AGN in the VCV c atalog is NGC 404 at D ~ 3 Mpc in the constellation 



Andromeda (IKarachentsev et al.ll2004l ). So the model AGNs with the distance from each 



observer D < D min = 3 Mpc were excluded. 

We checked the angular distance Q between a given reference object to its nearest 
neighboring object. For a set of 442 objects (the number of the AGNs with z < 0.018 in 
the VCV catalog), if they are distributed isotropically over the sky of Att radian, the average 
value of Q would be (Qiso) ~ H°- With the 442 AGNs from the VCV catalog, on the other 
hand, (Qvcv) = 3.55°. The fact that (Qvcv) < (Qiso) means that the distribution of the 
AGNs from the VCV catalog is not isotropic, but highly clustered, following the matter 
distribution in the LSS of the universe. We note that (Qvcv) = 3.55° is similar to the 
angular window of 3.1° used in the Auger study. Clearly this agreement is not accidental, 
but rather consequential. For the sets of ~ 500 model AGNs in our simulations, the average 
angular distance is (Qagn) = 3.68° ±1.66°. The error was estimated with (Qagn) for ~ 1400 
observers. That fact that (Qagn) ~ (Qvcv) indicates that the spatial clustering of our model 
AGNs is on average comparable to that of the AGNs from the VCV catalog. This provides 
a justification for our selection criteria for model AGNs in the simulated universe. We note 
that Q is an intrinsic property of the distribution of the reference objects in the sky and has 
nothing to do with UHECRs. 



2.5. Sources of UHECRs 



Although AGN is one of viable candidates that would produce UHECRs (see, e.g., 
Nagano fc Watson! |2000| . for the list of viable candidates), there is no compelling reason that 
all the nearby AGNs are the sources of UHECRs. In this paper, we considered three models 
with different numbers of sources, N STC (see Table 1), to represent different subsets of AGNs. 
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1) Among AGNs, radio galaxies are consid ered to be the most promising sources of UHECRs 



^see, e.g., iBiermann &: Strittmatterl 119871 ). and there are 28 known radio galaxies within 
D < 75 Mpc. So in Model C, we considered on average 28 model AGNs located at 28 hottest 
host clusters (kT x > 0.8keV) within a sphere of radius 75 Mpc as true sources of UHECRs. 2) 
Based on the ratio of singlet to doublet events, on the other hand, Abraham et al.l (j2008al ) 
argued that the lower limit on the number of sources of UHECRs would be around 61. 
Following this claim, in Model B, we regarded on average 60 model AGNs located at 60 
hottest host clusters (kT x > 0.55keV) as true sources of UHECRs. 2) In Model A, we 
regarded all the model AGNs (reference objects) as true sources of UHECRs. 



2.6. Simulations of Propagation of UHE Protons 



At sources, UHE protons were injected with power-law energy spectrum; N^E^) oc 
EZ? for 6 x 10 19 eV < E^ < 10 21 eV, where 7 is the injection spectral index. We considered 



the two cases of 7 = 2.7 and 2.4. At each source, protons were randomly distributed over a 
sphere of radius O.S/i^ 1 Mpc, and then launched in random directions. 

We then followed the trajectories of UHE protons in our model universe with the IGMF, 
by numerically integrating the equations of motion; 



dr 
~dt 



v. 



dp 

~dt 



e (v x B) 



(2) 



where p is the momentum. During propagation UHE protons interact with the CMB, and 
the dominant processes for energy loss are the pion and pair productions. The energy loss 
was treated with the continuous-loss approximation (IBerezinsky et al.ll2006l ). The adiabatic 
loss due to the cosmic expansion was ignored. 



Table 1. Models of different numbers of sources 

















Model 


N 

1 v src 


host clusters 




6 


('S'sim) 


9 ■ 


A 


-500 


kT x > O.lkeV 


13.98 


7.01 


3.58 


2.80 


B 


60 


kT x > 0.55keV 


15.33 


8.80 


3.97 


3.19 


C 


28 


kT x > 0.8keV 


17.76 


10.45 


4.23 


3.43 



a Deflection angle 9 and separation angle S are defined in Sec- 
tion 3 
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We let UHE protons continue the journey, visiting several observers during flight, until 
the energy falls to 60 EeV. At observers, the events with E > 60 EeV were recorded and 
analyzed. 

3. Results 

3.1. Deflection Angle 

In Paper I, we considered the deflection angle, 9, between the arrival direction of UHECR 
events and the sky position of their sources (see Figure 2). Obviously this angle can be cal- 
culated only when the true sources are known, which is not the case in experiments. In 
the simulated universe, our model AGNs and observers are located in strongly magnetized 
regions. As illustrated in Figure 1, UHECRs first have to escape from magnetic halos sur- 
rounding sources, then travel through more or less void regions (path 1) or through filaments 
(path 2), and finally penetrate into magnetic halos around observers, to reach observers. So 
the degree of deflection depends not only on the magnetic field along trajectories but also on 
the fields at host clusters and groups of sources and observers. Since the gas temperature, 
the depth of gravitational potential well, and the magnetic field energy density are related 
as kTx oc <3> oc in our model, hotter clusters and groups would have stronger fields. So 
we expect that if sources and observers are located at hotter hosts, 9 would be larger on 
average. 

Figure 3 shows 9 versus Dg for UHE proton events recorded at observers in our simula- 
tions. Here, Dg denotes the distance to the actual sources of events. The top, middle, and 
bottom panels are for Models A, B, and C, respectively. Only the case of injection spectral 
index 7 = 2.7 is presented. The case of 7 = 2.4 is similar. Each dot represents one simulated 
event, and there are about 10 5 events in each Model. The upper circles connected with 
dotted lines represent the mean values of 9 in the distance bins of [Dg, Dg + ADg]. The mean 
deflection angle averaged over all the simulated events is (9) = 13.98°, 15.33°, and 17.76° for 
7 = 2.7 in Model A, Model B, and Model C, respectively. The lower circles connected with 
solid lines represent the median values of 9. The median value for all the simulated events 
is = 7.01°, 8.80°, and 10.45° for 7 = 2.7 in Model A, Model B, and Model C, respectively. 
The values of (9) and 9 for 7 = 2.4 are similar. The marks connected with vertical solid 
lines on the both sides of median values are the first and third quartiles in the distance bins, 
which provide a measure of the dispersion of 9. 

We note the following points. (1) With our IGMF, the mean deflection angle of UHE 
protons due to the IGMF is quite large, much larger than the angular window of 3.1° used 
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in the Auger correlation study. It is also much larger than the mean deflecti on angle that is 
expec ted to result from the Galactic magnetic field, which is a few degrees (ITakami fc Sato 



20081 ). (2) The mean deflection angle is largest in Model C and smallest in Model A. Recall 



that in Model C sources are located only at 28 hottest clusters, while in Model A all 500 
clusters include sources (see Table 1). The UHECR events from hotter clusters tend to 
experience more deflection, as noted above. So the mean of deflection angles in the model 
with hotter host clusters is larger. (3) The mean value of 9 has a minimum at Dg min ~ 20 — 30 
Mpc, which compares to a typical length of filaments. As we pointed out in Paper I, this is 
a consequence of the structured magnetic fields that are concentrated along filaments and at 
clusters. In an event with Dg < Dg tm[n , the source and observer are more likely to belong to 
the same filament, and so the particle is more likely to travel through strongly magnetized 
regions and suffer large deflection (see the path 2 in Figure 1). In the opposite regime, the 
source and observer are likely to belong to different filaments, so the particle would travel 
through void regions (see the path 1 in Figure 1). (4) For the events with Dg > Dg tmin , the 
mean and dispersion of 9 increase with Dg. Such trend is expected, since in the diffusive 
transport model of the pro pagation of UHECRs, th e deflection angle increases with distance 



as # rms oc \JTTg (see, e.g., iKotera fc Lemoind 120091 . and references therein). (5) There are 



more events from nearby sources than from distant sources, although all the sources inject 
the same number of UHECRs in our model. The smaller number of events for larger Dg 
should be mostly a consequence of energy loss due to the interaction with the CMB. 



3.2. Separation Angles 

In this study, we also consider the separation angle, S, between the arrival direction of 
UHECR events and the sky position of nearest reference objects (see Figure 2). The angle can 
be calculated with observation data, once a class of refere nce objects (e.g. AGNs, galaxies 



gamm a-ray bursts, and etc.) is specified. For example, iThe Pierre Auger Collaboration 



( 120071 ) took the AGNs within 75 Mpc in the VCV catalog as the reference objects for their 
correlation study. However, for a given UHECR event, the nearest AGN in the sky may not 
be the actual source; hence, the separation angle between a UHEC R event and its nearest 



AGN is not necessarily the same as the deflection angle of the event (1 Hillaa 1 2 9l ; iRyu et al. 



2009|). 



We obtained S for simulated events with our model reference objects (AGNs). Figure 
4 shows S versus Ds- Here, Ds denotes the distance to nearest AGNs. Again only the 
case of 7 = 2.7 is presented, and the case of 7 = 2.4 is similar (see Figure 5). The circles 
connected with solid line represent the mean values of S for the events with nearest AGNs 
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in the distance bins of [Ds,D$ + A£)g]. The mean separation angle averaged over all the 
simulated events is (S sim ) = 3.58°, 3.97°, and 4.23° for 7 = 2.7 in Models A, B, and C, 
respectively. 

We note the following points. (1) The mean separation angle is much smaller than the 
mean deflection angle, (S 8 i m ) ~ (1/4) (9) (see the next section for further discussion). (2) 
The mean separation angle is largest in Model C and smallest in Model A, although the 
difference of (S B i m ) among the models is less than that of (9). With larger deflection angles 
in Model C, there is a higher probability for a event to be found further away from the region 
where model AGNs are clustered, so S is on average larger as well. (3) Similarly as in the 9 
versus Dg distribution, the distribution of (Ssim) has the minimum at around Ds ~ 35 Mpc. 
This is again a signature of the filamentary structures of the LSS. (4) Contrary to D e , there 
are more events with larger Ds than with smaller Ds- It is simply because there are more 
AGNs with larger Ds- 



3.3. Comparison with the Auger Data 



We also obtained S for the 27 Auger events of highest energies, published in Abraham et al 



( 12008al ). with 442 nearby AGNs from the VC catalog. In Figure 5, we compare the S versus 
D s distribution for the Auger data with that of our simulations. The upper circles connected 
with dashed/dot-dashed lines represent the mean values of S sirn of simulated events as in Fig- 
ure 4, but this time both cases of 7 = 2.7 and 2.4 are presented. The lower circles connected 
with solid/dotted lines represent the median values of S S i m . The difference between the cases 
of 7 = 2.7 and 2.4 is indeed small. The median value of S B i m for all the simulated events 
is S sim = 2.80°, 3.19°, and 3.43° for 7 = 2.7 in Models A, B, and C, respectively. Asterisks 
denote the Auger events. The mean separation angle for the Auger events is (^Auger) = 3.23° 
for 26 events, excluding one event with large S (~ 27°), while (SAuger) = 4.13° for all the 
27 events. We note that (^Auger) ~ (Ssim), even though the mean deflection angle is much 
larger than the mean separation angle in our simulations, that is, (9) ^> (S 8 i m ). In all the 
models, about a half of the Auger events lie within the quartile marks: 15, 13, and 13 events 
for Models A, B, and C, respectively. 

With (9) ~ 15° in our simulations, one might naively expect that such large deflection 
would erase the anisotropy in the arrival direction and the correlation between UHECR 
events and AGNs (or the LSS of the universe). However, we argue that the large deflection 
does not necessarily lead to the general isotropy of UHECR arrival direction, if the agent of 
deflection, the IGMF, traces the local LSS. Suppose that UHECRs are ejected from sources 
inside the Local Supercluster. Some of them will fly along the Supergalactic plane and arrive 
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at the Earth; their trajectories would be deflected by the magnetic field between sources and 
us, but the arrival directions still point toward the Supergalactic plane. Others may be 
deflected into void regions, and then they will have less chance to get reflected back to the 
direction toward us due to lack of the turbulent IGMF there. In a simplified picture, we may 
regard the irregularities in the IGMF as the 'scat ters' of UHECRs; then t he last scattering 



point will be the arrival direction of UHECRs (see lKotera fc Lemoindl2009l . for a description 
of deflection of UHECRs based on this picture). As a result, even with large deflection, we 
still see more UHECRs from the LSS of clusters, groups, and filaments, and fewer UHECRs 
from void regions where both sources and scatters are underpopulated. Consequently, the 
anisotropy in the arrival direction of UHECRs can be maintained and the arrival direction 
still follows the LSS of the universe. 

Below the GZK energy, the proton horizon reaches out to a few Gpc, so the source 
distribution should look more or less isotropic and the arrival directions should not show a 
correlation with nearby AGNs. Thus, we do not expect to see anisotropy and correlation for 
UHECRs with such energy. 

In Section 2.4, we showed that the degree of clustering of our model AGNs is similar to 
that of AGNs from the VCV catalog; the mean of the angular distance Q between a given 
AGN to its nearest neighboring AGN is similar, (Qagn) ~ (Qvcv)- In both cases, AGNs 
follow the matter distribution in the LSS, highly structured and clustered. We point that 
if along with the reference objects, the CR sources and the IGMF also follow the matter 
distribution, with (9) ^> (S), (S) ~ (Q) is expected. The result that (SAuger) ~ (Ssim) 
~ (Qagn) ~ (Qvcv) is indeed consistent with such expectation. This means, however, that 
the statistics of S reflect mainly on the distribution of reference objects, rather than the 
deflection angle. 

To further compare the Auger data with our simulations, we plot the cumulative fraction 
of events, F(< log S), versus logS for the simulated events (lines) and the Auger events 
(open circles) in Figure 6. The solid and dotted lines are for the cases of 7 = 2.7 and 
2.4, respectively, and the difference between the two cases is again small. The Kolmogorov- 
Smirnov (K-S) test yields the maximum difference of D = 0.17, 0.23, and 0.26 between the 
Auger data and the simulation data (7 = 2.7) in Models A, B, and C, respectively; the 
significance level of the null hypothesis that the two distributions are statistically identical 
is P ~ 0.37, 0.09, and 0.04 for Models A, B, and C, respectively. So the null hypothesis 
that the two distributions for our simulated events and the Auger events are statistically 
identical cannot be rejected, especially for Model A. This would be a justification for our 
models of the IGMF, sources of UHECRs, and reference objects. Also we see that Model A 
with more sources is preferred over Models B and C with fewer sources. But this does not 
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necessary mean that all the AGNs would be the actual sources of UHECRs. We note that 
the number of the Auger events used, 27, is still limited. In addition, we consider only UHE 
protons in this paper. Hence, before we argue the above statements for sure, we will need 
more observational events and need to know the composition of UHECRs (see Summary and 
Discussion for further discussion on composition). 

3.4. Probability of Finding True Sources 

With (9) ^> (S S i m ), there is a good chance that the AGNs found closest to the direction 
of UHECRs are not the actual sources of UHECRs. To illustrate this point, we first show the 
distribution of D$ versus Dg in Figure 7. For some events, the closest AGNs are the actual 
sources. They are represented by the diagonal line. Around the diagonal line, a noticeable 
fraction of events are found. Those are the events for which the closest AGNs are found 
around the true sources; both sources and close-by AGNs are clustered as a part of the LSS 
of the universe. For the events away from the diagonal line, it is more likely that D s > D e . 
It is because there are more AGNs with larger Ds; away from true sources, observed events 
are more likely to pick up closest AGNs with larger D S - 

To quantify the consequence of (9) ^> (S siin ), we calculated the fraction of true identifi- 
cation, /, as the ratio of the number of events for which nearest AGNs are their true sources 
to the total number of simulated events. This is a measure of the probability to find the 
true sources of UHECRs, when nearest candidates are blindly chosen (which is the best we 
can do with observed data). In Figure 8, we show the fraction as a function of separation 
angle, S. The fraction is largest in Model A with largest N STC , and smallest in Model C with 
smallest N src . At S ~ 2° the fraction is about 50 % for Model A, close to 40 % for Model B, 
and a little above 30 % for Model C, but only 20 - 30 % at S = 3° - 4°. As the separation 
angle increases, the fraction decreases gradually to ~ 10 %, indicating lower probabilities to 
find true sources at larger separation angles. On average, we should expect that in less than 
1 out of 3 events, the true sources of UHECRs can be identified, if our model for the IGMF 
is valid and UHECRs are protons. 

4. Summary and Discussion 

In the search for the nature and origin of UHECRs, understanding the propagation of 
charged particles through the magnetized LSS of the universe is important. At present, the 
details of the IGMF are still uncertain, mainly due to limited available information from 
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observation. Here, we adopted a realistic model universe that was described by simulations 
of cosmological structure formation; our simulated universe represents the LSS, which is 
dominated by the cosmic web of filaments interconnecting clusters and groups. The distri- 
bution of the IGMF in the LSS of t he universe was obtained with a physically motivated 



model based on turbulence dynamo (IRyu et al.ll2008l ). 



To investigate the effects of the IGMF on the arrival direction of UHECRs, we further 
adopted the following models. Virtual observers of about 1400 were placed at groups of 
galaxies, which represent statistically the Local Group in the simulated model universe. 
Then, we set up a set of about 500 AGN-like "reference objects" within 75 Mpc from each 
observer, at clusters of galaxies (deep gravitational potential wells) along the LSS. They 
represent a class of astronomical objects with which we performed a correlation analysis for 
simulated UHECR events. We considered three models, in which subsets of the reference 
objects were selected as AGN-like sources of UHECRs (see Table 1). UHE protons of E > 60 
EeV with power-low energy spectrum were injected at those sources, and the trajectories of 
UHE protons in the magnetized cosmic web were followed. At observer locations, the events 
with E > 60 EeV from sources within a sphere of radius 75 Mpc were recorded and analyzed. 

To characterize the clustering of the reference objects, we calculated the angular dis- 
tance, Q, from a given reference object to its nearest neighbor. The mean value for our 
model AGNs in the simulated universe is (Qagn) = 3.68° ± 1.66°, while that for 442 AGNs 
from the VCV catalog is (Qycv) — 3.55°. This demonstrates that the two samples have a 
similar degree of clustering and are highly structured (e.g. (Qi BO ) ~ 11° for the isotropic 
distribution). 

With our model IGMF, the deflection angle, 9, between the arrival direction of UHE 
protons and the sky position of their actual sources, is quite large with the mean {9) ~ 14 — 
17.5° and the median 9 ~ 7 — 10°, depending on models with different numbers of sources (see 
Table 1). On the other hand, the separation angle between the arrival direction and the sky 
position of nearest reference objects is substantially smaller with the mean (S S i m ) ~ 3.5 — 4° 
and the median S sim = 2.8 — 3.5°. That is, we found that while (9) ~ 4( < S , siTT1 ), (^Jis 



similar to (Qagn)- F° r the Auger events of highest energies in Abraham et al.l (l2008al ). 
with 442 nearby AGNs from the VCV catalog as the reference objects, the mean separation 
angle is (^Auger) = 3.23° for the 26 events, excluding one event with large S (~ 27°), while 
(^Auger) = 4.13° for all the 27 events. Hence, (SAuger) ~ (Qycv) ~ (S s i m ) ~ (Qagn)- This 
implies that the separation angle from the Auger data would be determined primarily by the 
distribution of reference objects (AGNs), and may not represent the true deflection angle. 

We further tested whether the distributions of separation angle, S, for our simulated 
events and for the Auger events are statistically comparable to each other. According the 
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Kolmogorov-Smirnov test for the cumulative fraction of events, F(< log S), versus logS 1 , 
the significance level of the null hypothesis that the two distributions are drawn from the 
identical population is as large as P ~ 0.37 for Model A (see Table 1). Thus, we argued that 
our simulation data, especially in Model A, are in a fair agreement with the Auger data. 
This test also showed that the model with more sources (Model A) is preferred over the 
models with fewer sources (Models B and C). 

The fact that (9) ^> (S s i m ) implies that the AGNs found closest to the direction of 
UHECRs may not be the true sources of UHECRs. We estimated the probability of finding 
the true sources of UHECRs, when nearest reference objects are blindly chosen: f(S) is the 
ratio of the number of true source identifications to the total number of simulated events. 
This probability is ~ 50 — 30 % at S ~ 2°, but decreases to ~ 10 % at larger separation angle. 
On average, in less than 1 out of 3 events, the true sources of UHECRs can be identified in 
our simulations, when nearest reference objects are chosen. 

The distribution of 9 versus Dq shows a bimodal pattern in which 9 is on average larger 
either for nearby sources (for Dq < 15 Mpc) or for distant sources (for Dg > 30 Mpc) with 
the minimum at the intermediate distance of -D^min ~ 20 — 30 Mpc. The distribution of 
S versus D s shows a similar, but weaker sign of the bimodal pattern. This behavior is a 
characteristic signature of the magnetized cosmic web of the universe, where filaments are 
the most dominant structure. When a large number of super-GZK events are accumulated, 
we may find the signature of the cosmic web of filaments in the S versus D s distribution. 

Finally, we address the limitations of our work. (1) We worked in a simulated universe 
with specific models for the elements such as the IGMF, observers, sources, and reference 
objects, but not in the real universe. So we could make only statistical statements. (2) 
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argue that our model for the IGMF is most plausible, since it is a physically motivated mode l 



based on turbulence dynamo without involving an arbitrary normalization (IRyu et al.ll2008l ). 
Nevertheless, our IGMF model should be confirmed further by observation. (3) The sources 
of UHECRs may not be objects like AGNs, but could be objects extinguished a while ago 



such as gamma-ray bursts (see, e.g., Vietri fl995l; Waxman Il995l). or sour ces spread over 



space like cosmological shocks (see, e.g., Kang et. al. 199(k Kang et al. 1997 ). The injection 
energy spectrum of power-law with cut-off at an arbitrary maximum energy (see Section 2.6) 
would be unrealistic. The IGMF in the Local Group (see Paper I), although currently little 
is known, might be strong enough to substantially deflect the trajectories of UHECRs. All of 
those will have effects on the quantitative results, which should be investigated further. (4) 
Recently, the Auger collaboration disclosed the analysis, which suggests a substantial fraction 
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of highest energy UHECRs might be iron nuclei (jUnger et all 120071 ; IWahberg et all 120091 ) . 



This is in contradiction with the an alysis of the HiRes data, wh ich indicates highest energy 



UHECRs would be mostly protons (jSokolsky &: Thomson! 120071 ). The issue of composition 



still needs to be settled down among experiments. Iron nuclei, on the way from sources to 
us, suffer much larger deflection than protons. Hence, if a substantial fraction of UHECRs is 
iron, some of our findings will change, a question which should be investigated in the future. 
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Fig. 1. — Distribution of the IGMF in a two-dimensional slice of (143 Mpc) 2 in the simulated 
universe. Locations of virtual observers (circles) and model AGNs (stars) are schematically 
marked at clusters and groups of galaxies. Paths of UHECRs from sources to observers are 
also schematically drawn. 
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Fig. 2. — Geometry of the deflection angle, 9, the separation angle, S, the distance to the 
true source, D e , and the distance to the nearest object, D s . The path of the UHECR event 
from the source to observer is schematically drawn. 
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Fig. 3. — Distribution of deflection angles (8) as a function of distance to true sources {Do). 
Dots represent UHE proton events recorded in our simulations. Circles connected with 
dotted lines (upper) and solid lines (lower) show the mean and median values, respectively. 
Vertical lines connect the marks of first and third quartiles in given D 9 bins. Top, middle, 
and bottom panels are for Model A, Model B, and Model C, respectively. The cases of 
7 = 2.7 are shown. 
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Fig. 4. — Distribution of separation angles between the directions of UHE protons and 
nearest model AGNs (S) as a function of distance to nearest model AGNs (-D5). Dots 
represent UHE proton events recorded in our simulations. Circles connected with solid lines 
show the mean values. Top, middle, and bottom panels are for Model A, Model B, and 
Model C, respectively. The cases of iniection spectral index 7 = 2.7 are shown. 
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Fig. 5. — Mean (upper circles) and median (lower circles) values of S as a function of Ds for 
UHE proton events in our simulations. Vertical lines connect the marks of first and third 
quartiles in given Ds bins. Circles connected with dashed lines (upper) and solid lines (lower) 
are for 7 = 2.7, and those connected with dot-dashed lines (upper) and dotted lines (lower) 
are for 7 = 2.4. The mean S for 7 = 2.7 are the same as those in Figure 2. The median 
and quartiles for 7 = 2.4 are horizontally shifted for b etter visibility. Asterisk s denote S for 
the 27 Auger events of highest energies, published in lAbraham et al.l (j2008al ). with nearby 
AGNs from the VCV catalog. Top, middle, and bottom panels are for Model A, Model B, 
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Fig. 6. — Cumulative fraction of UHECR events with separation angle smaller than S. 
Solid and dotted lines denote the results calculated with the UHE proton events in our 
simulations for 7 = 2.7 and 2.4, respectively. Top, middle, and bottom panels are for Model 
A, Mo del B, and Model C, r espectively. Open circles denote the result for the 27 Auger 
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Fig. 7. — Distance to nearest model AGNs, D s , versus distance to true sources, D e , for the 
UHECR events in Model A. The case of 7 = 2.7 is shown. Color codes the number of events 
in log 10 scale in bins of AD e x AD S - The maximum distance is 75 Mpc for both D s and D g . 
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Fig. 8. — Fraction (/) of UHE proton events recorded in our simulations, for which their 
true sources are identified as closest AGNs in the sky, as a function of S. Dashed, dotted, 
and solid lines are for Model A, Model B, and Model C, respectively. Heavy and light lines 
denote the fractions for 7 = 2.7 and 2.4, respectively. 



